HIGH-SPEED LINE-BY-LINE CALCULATION METHOD AND 
PROGRAM BY DIVISION OF VOIGT PROFILE 



TECHNICAL FIELD 

[0001] The present invention relates to high-speed line-by-line calculation of the absorption 
coefficient spectra of gas molecules necessary for atmospheric observations. 

BACKGROUND ART 

[0002] While the Earth's atmosphere receives electromagnetic radiation from the sun, it is 
itself releasing electromagnetic radiation mainly in the infrared region into space. This 
electromagnetic radiation interacts such as by absorption, radiation and scattering with gas 
molecules, clouds, aerosols and the like present in the atmosphere. By measuring changes in 
the spectral distribution of electromagnetic radiation caused by such interactions with the 
atmosphere (theoretical calculations of such distributions are referred to as radiative transfer 
calculations) using sensors on land or in airplanes or satellites, it is possible to study the air 
temperature, air pressure and concentrations of atmospheric trace gas constituents. 
[0003] While a techruque known as line-by-line calculation is used for radiative transfer 
calculations of the atmosphere, line-by-line calculations involve a enormous number of 
calculations. For this reason, the calculation time becomes extremely long, making it 
practically difficult to perform calculations over a wide spectral range at high speeds using 
personal computers or the like. Additionally, there are methods that make use of data 
obtained in line-by-line calculations beforehand involving table lookup routines as methods 
for avoiding the problem discussed above, but such methods have the drawback of requiring 



large data files. Additionally, other algorithms not requiring prior calculation have fared no 
better than about 20 times the speed of line-by-line calculations. 

Non-Patent Publicaiton 1: Akihiro Uchiyama, "Line-by-Line Computation of the 
Atmospheric Absorption Spectrum Using the Decomposed VOIGT Line Shape", /. Quant, 
Spectrosc. Radiat. Transfer, Vol. 47, No. 6, pp. 521-532, 1992. 

DISLCOSURE OF THE INVENTION 
Problems to be Solved by the Invention 

[0004] Due to the problems mentioned above, calculation methods capable of line-by-line 
calculations at high speed have been sought, and the present invention offers such a 
Ccdculation method and program. 

Means for Solving the Problems 

[0005] The present invention relations to a Voigt function approximation calculating 
program used for line-by-line calculations, said program performing: 

(1) a step of dividing a domain of a Voigt function into a first range around the peak of 
the Voigt function and a skirt portion not contained in the first range, replacing the first 
range with a cubic function, calculating the values and derivatives of said cubic function 
and the Voigt function in the skirt portion for each of first predetermined intervals, and 
connecting said cubic function and said Voigt function at points of connection thereof 
using the values and derivatives of both functions; 

(2) a step of adding together the results of step (1) for a plurality of absorption lines; 

(3) a step of calculating the function values and derivatives for the results of step (2) by 
interpolation over intervals smaller than said first predetermined intervals; 
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(4) a step of dividing said first range into a second range near the peak and a skirt 
portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
cubic function, and calculating values and derivatives of said cubic function and said 
"function representing the difference between the Voigt function and said cubic function" 
in the skirt portion for each of second predetermined intervals; 

(5) a step of connecting Scdd cubic function and said "function representing the 
difference between the Voigt function and said cubic function" at points of connection 
thereof using the values and derivatives of both functions; 

(6) a step of adding the results of steps (4) and (5) to the results of step (3) for a plurality 
of absorption lines; 

(7) a step of calculating the function values and derivatives for the results of step (6) by 
interpolation over intervals smaller than said second predetermined intervals; and 

(8) a step of adding the values of the "function representing the difference between the 
Voigt function and said cubic function" to the results of step (7) for a plurality of 
absorption lines in said second range. 

A method using the above steps allows a Voigt function to be calculated at high 
speed using a cubic function. 

[0006] In addition thereto, the steps (5) through (8) may be repeated one or more times imtil 
a third predetermined interval is reached. As a result, it is possible to obtain detailed results 
with small calculation intervals. The first through third predetermined intervals can be 
determined as described below. 

[0007] The first predetermined interval for the widest sub-function is j^^dv. Here, ; is a 
single-digit natural number, dv is the increment in wave number, and kmax is the largest 
natural number satisfying the relationship j^^^^^pdv < Vmax, However, Vmax represents the 
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maximum calculation range from the center of the absorption line, and p is a natural number (p 
= 1, 2, 3) for controlling the calculation precision 

The most specific third predetermined interval is j^*"dv. Here, ; is a single-digit 
natural number of about 2-6, dv is the increment of the wave number, and kmin is the 
maximum non-negative decimal fraction (0 when non-existent) satisfying the relationship 
jkminp^y <a(a being approximately y/4). However, y is an approximate value of the full-v^ddth 
at half-max of the absorption line, and p is a natural number (p = 1, 2, 3) for controlling the 
calculation precision. 

A method in accordance with any one of claims 1-4, wherein the second 
predetermined intervals are determined using the following equation. 

For the sub-function with the (k - kmin + l)-th smallest vddth, the predetermined 
interval is /rfv. Here, / is a single-digit natural number, dv is the increase in wave number and 
k is such that kmin <k< kmax, 

[0008] For a quartered interpolation, the interpolation is calculated with; set to 4, using the 
function values yo, yi and function derivative values yo', yi' at xo, xi in the interpolation interval 
{xo, xi), with the below-given Equation (1) as a function value interpolation equation. Equation 
(2) as a function derivative value interpolation equation and £ as a non-negative decimal 
fraction. 
[Eq. 1] 
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[0009] For an interpolation divided into five parts, the interpolation is calculated with ; set 
to 5, using the function values yo, yi and function derivative values yo', yi' at xo, xi in the 
interpolation interval (xo, xi), with the below-given Equation (3) as a function value 
interpolation equation. Equation (4) as a function derivative value interpolation equation and e 
as a non-negative decimal fraction. 
[Eq.31 
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[0010] Additionally, the invention offers a method for achieving high speeds, when 
£issuming the Voigt function to be K(x, y) and ttie difference from ttie Voigt profile of the 
absorption line to be K{x, y) +J[x), by replacing 
[Eq.5] 

K{x,y)=AKix,y)+Bf{x) 
and 
[Eq. 6] 
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dK{x,y) 
dx 

[0011] Furthermore, the invention offers a method for achieving high speeds, when 
assuming the difference from the Voigt profile to be K(x, y)f(x), by replacing 
[Eq. 7] 

K{x,y) = K{x,y)f{x) 
and 
[Eq. 8] 

^E^ = ^^fi:c)^K{x,y)^ 
dx dx ■'^ ' ^ dx 

[0012] Additionally, in line-mixing correction: 

[Eq.9] 

K{x,y) = AK(x,y)+ BL{x,y) 
and 



[Eq. 10] 
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Here, L{x, y) is the imaginary component of the function w{z) (the real part is the Voigt 
function), where complex number z = x + \y, defined by the following equation. 
[Eq. 11] 

w(z) = ^ £ ^^^ti!) dt = exp(- )erfc{- iz) = K{x, y) + iL{x, y) 
{erfc{z) is a complex complementary error function). 

[0013] Furthermore, the invention offers a device for performing the method of the present 
invention by loading the method described above as a program. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



[0014] [Fig. 1] Fig. 1 is a schematic view of a Voigt function. 

[Fig. 2] Fig. 2 is a diagram showing how a Voigt function is divided into 
sub-functions according to the method of the present invention. 

[Fig. 3] Fig. 3 is a diagram showing the error between a Voigt function and a cubic 

function. 

Explanation of Reference Numbers 
[0015] 1 Voigt function 

2 first range 

3 skirt portion 

4 function contained in first range 

5 portion not contained in first range 



BEST MODE FOR CARRYING OUT THE INVENTION 

[0016] An embodiment of the present invention shall be described below. Fig. 1 is a 
diagram showing the shape of a Voigt function 1. According to the present invention, this 
Voigt function 1 is first divided into sub-functions. These sub-functions are in a first range 2 
around the peak and a skirt portion 3 outside the first range 2. This first range 1 is replaced 
by a cubic function, and the portion not contained in the first range, in other words, the skirt 
portion, is calculated with the Voigt function. 

[0017] The details of the calculation method are shown in Fig. 2. As shown in Fig. 2(a), the 
function 4 ,which is the portion of the Voigt function 1 contained in the first range 2, is 
replaced by a cubic function. At this time, the cubic function is determined so as to smooth 
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the connections by matching the derivatives of the functions at the points of connection 
between the cubic function and the Voigt function. Additionally, the portion 5 not contained 
in the first range is calculated by the Voigt function. The calculation is performed over a first 
predetermined interval, indicated by j^^^dv where dv is the change in wave number. Here, ;* is 
a single-digit natural number, and kmax is the largest natural number satisfying the 
relationship j^^^pdv < Vmax. Vmax may, for example, be 25 cm Additionally, ; is normally 
a number of about 2-6, and is preferably 4. Also, Vmax represents the calculation range from 
the center of the absorption line, and p is a natural number for controlling the calculation 
precision, p preferably having a value of 1, 2 or 3. 

[0018] Results calculated in this way are added together for a plurality of absorption lines. 
At this stage, the result is an approximate value calculated in the first predetermined interval. 
[0019] Next, unevaluated points are calculated over an interval smaller than the first 
predetermined interval. That is, the unevaluated points are interpolated over intervals 
obtained by quartering the first predetermined interval. The calculation of the point of 
division is performed by the following interpolation formula. 

The function values yo, yi and function derivative values yo', yi at xo, xi in the 
interpolation interval (xo, xi) can be determined by the following two equations: 
[Eq. 12] 
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Here, £ is a non-negative decimal fraction. With regard to e, even if the Voigt function and 
cubic function match at the interpolation points xo, xi as shown in Fig 3, there will be an error 
in the function values due to differences in the curvature of the curves in the region between 
the interpolation points. While adequate precision can be obtained even if £ = 0, the error can 
be further reduced by about half by changing the value of e to adjust the curvature of the cubic 
function. From these two equations it is possible to determine the function values ya, yb, yc 
and the derivative values ya, yb', y/ at points where the interpolation interval (xo, xi) of the 
aforementioned first predetermined interval has been quartered. Here, xo, xi are wave 
numbers in the present invention. The aforementioned adjustment of e is performed by 
comparing the calculation results when changing £ with the calculation results for a normal 
line-by-line method for a small spectral range which is only a part of the spectral range to be 
calculated. 

[0020] Additionally, the present invention proposes an interpolation in which the first 
predetermined interval is divided into five parts. The interpolation equations for dividing 
into five parts are as follows: 
[Eq. 14] 
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The symbols used in the equations are the same as in the case of quartering described above. 
Additionally, those skilled in the art would recognize that yd, y/ clearly represent a function 
value and a derivative value. Here, the interpolation equation described above can be used 
not only for the first predetermined interval, but also for other predetermined intervals, and 
this should also be obvious to those skilled in the art. 

[0021] The interpolation points are calculated as described above. Next, as shown in Fig. 

2(b), further calculations are performed by interpolation of the region aroimd the peak in the 

first range. As described above, the area around the peak in the first range is called the 

second range 6, and divided from a skirt portion not contained in the second range 6. Then 

the second range is represented by a second cubic function, and the skirt portion not contained 

in the second range is represented by a "function representing the difference between the 

Voigt function and the first cubic function" 7. This second cubic function represents the 

difference between the new cubic function and the first cubic function. Here, this "function 

representing the difference between the Voigt function and the first cubic function" shall 

henceforth be referred to simply as the "difference function". As described above, the second 

cubic function and the "difference function" 7 are again joined under the condition that at the 

points of connection, the function values and derivative values are matched. 

[0022] After connecting the functions in this way, the function value and derivative value 

are calculated in the second predetermined interval for the second range. This second 

predetermined interval is smaller than the first predetermined interval and is calculated by the 

equation given below. 

For the sub-function with the (k - kmin + l)-th smallest width, the predetermined 
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interval is j^dv. Here, ; is a single-digit natural number, dv is the increase in wave number and 
k is such that kmin <k< kmax. 

Additionally, these calculations are performed for all absorption lines. Furthermore, 
interpolation is performed between the calculated points. The equations (1) and (2) or (3) and 
(4) are used for the interpolation. 

[0023] Here, the above method for calculating a "difference function" for a skirt portion not 
contained in the second range and a cubic function for the second range that connects 
smoothly therewith for a plurality of absorption lines has been explained, but it is also possible 
to use a method of calculating the Voigt function of the skirt portion not contained in the 
second range and a cubic function connecting smoothly therewith for a plurality of absorption 
lines, then subtracting the first cubic function from the entire first range for a plurality of 
absorption lines. 

[0024] With regard to this second range, the above calculations may be repeated until a 
predetermined interval is reached. That is, the second range may be further divided into a 
peak area and a skirt portion which are respectively calculated by a cubic function and 
difference function, then connected, and the calculation performed with k in the equation for 
calculating the second predetermined interval incremented by 1. This calculation is 
performed for all absorption lines, and these are added to the results calculated in Fig. 2(b). 
The steps described above are repeated until the predetermined interval is reached. 
[0025] Further, a calculation is performed after replacing the peak area with a cubic function 
and making the skirt portion a "difference function", over a third predetermined interval 
smaller than the second interval. This third predetermined interval is represented by the 
equation given below. 

The third predetermined interval is;^'«'"dv. Here, / is a single-digit natural number, 
dv is the increment of the wave number, and kmin is the maximum non-negative decimal 
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fraction satisfying the relationship j^*npdv < a (if non-existent, it is set to 0). However, a has a 
value of about y/4 (y being an estimate of the full width at half max of the absorption line), and 
p is a natural number (p = l,2, 3) controlling the calculation precision. 

[0026] By following the above steps, it is possible to calculate a Voigt function precisely and 
at high speed. Here, it is known that there is a slight difference from the Voigt function of 
carbon dioxide and the like, and a method for correcting for this difference is known. The 
method of the present invention can also be applied to such corrections. These corrections 
shall be explained next. 

[0027] General profile correction shall be explained first. 

The Voigt function is represented by K{x, y). When the difference from the Voigt 
profile is represented by the shape of K{x, y) + f{x), the high-speed techruque can be applied by 
setting: 
[Eq. 16] 

as the function to be substituted, and replacing it with: 
[Eq. 17] 

k{x,y)=AK{x,y)+Bf{x) 
[Eq. 18] 

dK{x,y) _ ^ dK{x,y) ^ ^dfjx) 
dx dx dx 

Similarly, if the difference from the Voigt profile is represented by the shape of K{x, y)f{x), then 

the high-speed technique can be performed by replacing it with 

[Eq. 19] 

K{x,y)^K{x,y)f{x) 
[Eq. 20] 
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OX OX dx 

[0028] Explaining the above correction method with a specific example, in a sub-Lorentzian 
correction: 
[Eq.21] 

[Eq. 22] 

dKix,y) 
dx 

will respectively become: 
[Eq. 23] 

K{x,y) = K{x,y)Acxp(-B\:^) 
[Eq. 24] 

Here, A and B are correction coefficients. Additionally, those skilled in the art will recognize 
that sgn(a:) is the sign function. 

[0029] For line-mixing correction, the following relationships arise: 
[Eq.25] 

^^^ = 2yL{x,y)-xKix,y)] 



[Eq. 26] 
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dx 

Therefore, 
[Eq. 27] 
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K{x,y) 
[Eq. 28] 

dK{x,y) 
dx 

will respectively be 
[Eq. 29] 

K{x,y)= AK{x,y) + BL{x,y) 
[Eq. 30] 

{Ax + By)K{x,y)- {Ay - Bx)L{x,y)-^ 

Here, L(x, y) is the imaginary component of the hmction w{z), where complex number z = x + iy, 

defined by the following equation. The real part is the Voigt function. 

[Eq.31] 

vi<z) = ^ L^^^^^dt = exp(- )erfc{- iz) = K{x,y)^ iL{x,y) 
erfc{z) is a complex complementary error function. 

[0030] According to the calculation methods of the present invention described above, it is 
possible to perform line-by-line calculations using Voigt functions with good precision and at 
high speed. By further loading these methods as computer programs, it is possible to result 
in specific devices. Such a program may, for example, be a stand-alone type program 
operating on a personal computer. It may also be provided in the form of a plug-in for other 
software. 
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